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1 Abstract 

We present results of a high statistics calculation of hadron masses and meson decay constants 
in the quenched approximation to lattice QCD with Wilson quarks at (3 — 5.85 and 6.0 on 
24 3 x 54 lattices. We analyze the data paying attention in particular to the systematic errors 
due to the choice of fitting range and due to the contamination from excited states. We find 
that the systematic errors for the hadron masses with quarks lighter than the strange quark 
amount to 1 — 2 times the statistical errors. When the lattice scale is fixed from the p 
meson mass, the masses of the Q~ baryon and the meson at two /3's agree with experiment 
within about one standard deviation. On the other hand, the central value of the nucleon 
mass at (3 = 6.0 (5.85) is larger than its experimental value by about 15% (20%) and that of 
the A mass by about 15% (4%): Even when the systematic errors are included, the baryon 
masses at (3 = 6.0 do not agree with experiment. Vector meson decay constants at two values 
of f3 agree well with each other and are consistent with experiment for a wide range of the 
quark mass, when we use current renormalization constants determined nonperturbatively 
by numerical simulations. The pion decay constant agrees with experiment albeit with large 
errors. Results for the masses of excited states of the p meson and the nucleon are also 
presented. 



1 Introduction 



Although there have been many efforts to calculate hadron masses in lattice QCD by nu- 
merical simulations, it has turned out that derivation of convincing results is much harder 
than thought at the beginning, even in the quenched approximation. For example, before 
1988, there was large discrepancy among the results for the mass ratio m^jm p obtained for 
(3 = Q/g 2 = 5.7 — 6.0 and in the quark mass region corresponding to m^jm p > 0.5. The 
discrepancy was caused by systematic errors due to contamination from excited states |T|, |2j 
and effects of finite lattice spacing |3j and finite lattice volume. Recent high statistics sim- 
ulations employ lattices with large temporal extent ||, [5|, |B| and/or extended quark sources 
H ||, [7], H |9], [10], [Tl| to reduce fluctuations as well as the contamination from excited states. 
However, a long plateau in an effective mass is rarely seen and data of effective masses 
frequently show large fluctuations at large time separations. The uncertainty in the choice 
of fitting range is therefore another source of systematic errors. In order to obtain reliable 
values for the spectrum, it is essential to make a quantitative study on these systematic 
errors. 

In this paper we report results of a high statistics calculation of the quenched QCD spec- 
trum with the Wilson quark action at (3 = 5.85 and 6.0 on 24 3 x 54 lattices. Our major 
objection is to calculate light hadron masses as well as meson decay constants paying atten- 
tion in particular to the systematic errors due to the choice of fitting range and due to the 
contamination from excited states. In order to estimate the magnitude of these systematic 
errors, we perform correlated one-mass fits to hadron propagators systematically varying 
fitting ranges || |i~2| . Assuming the ground state dominance at large time separations, we 
estimate systematic errors in hadron masses which cannot be properly taken into account by 
the standard least mean square fit when the fitting range is fixed. It is shown that, for the 
hadron masses with quarks lighter than the strange quark, the systematic errors amount to 
1 — 2 times the statistical errors. We then perform correlated two-mass fits, again varying 
fitting ranges. We find that the ground state mass is consistent with that obtained from the 
one-mass fit within the statistical and systematic errors. Finally, we extrapolate the results 
of hadron masses at finite quark mass to the chiral limit, taking account of systematic errors 
both due to the choice of extrapolation function and due to the fitting range. We also study 
meson decay constants in a similar way. 

We use the point source in this study. Historically there was a report that numerical 
results for hadron masses appear to depend on the type of the source adopted [[L3|] , although 
it has afterward reported in some works that masses are independent within the statistical 
errors || ||. Note in this connection that there is no proof that the value of a hadron mass 
is independent on the type of sources in the case of the quenched approximation due to the 
lack of the transfer matrix and that there is the so-called Gribov problem for gauge fixing 
which is necessary for almost all smeared sources. Under these circumstances it may be 
worthwhile to present the details of the results and the analyses with the point source as a 
reference. The method of analyses of the systematic errors in this work can be applied to 
the cases of smeared sources too. 

Numerical simulations are performed with the QCDPAX ||14|| , a MIMD parallel computer 
constructed at the University of Tsukuba. For the calculations performed in this work, we 
use 24 x 18 processing units interconnected in a toroidal two-dimensional mesh with a peak 
speed of 12.4 GFLOPS. (The maximum number of nodes is 24 x 20 with a peak speed of 14.0 
GFLOPS.) The sustained speed for the Wilson quark matrix multiplication is approximately 
5 GFLOPS. The calculations described here took about six months on the QCDPAX. 
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We start by giving in Sec. |2| some details about our numerical simulations. Then we 
derive hadron masses at finite quark mass in Sec. |3] and perform two-mass fits to estimate 
the masses of excited states of the p meson and the nucleon in Sec. |[ We extrapolate the 
results to the chiral limit in Sec. ^. Sec. [] is devoted to the evaluation of meson decay 
constants. In Sec. 0, we give conclusions and discussion on the results. 

2 Numerical Calculation 

We use the standard one-plaquette gauge action 

S g = ^BBTt{U P ) W 
9 p 

and the Wilson quark action 

S q = -J2Hn)D(K,n,m)^(m), (2) 

n,m 

D(K } n, m) = <5„, m - K ii 1 ~ T/J^A+m,™ + ( J + ^) u L,^ S m+ix,n}, (3) 

where g is the bare coupling constant and K is the hopping parameter. 

Simulations are done on 24 3 x 54 lattices at = 6/g 2 = 5.85 and 6.0 for the five 
values of hopping parameter listed in table III. The mass ratio m 7T /m p takes a value from 
0.97 to 0.52 and roughly agrees with each other at two /3's for the five cases of hopping 
parameter. We choose the values of the third largest hopping parameter in such a way that 
they approximately correspond to the strange quark. 

We generate 100 (200) configurations with periodic boundary conditions at (5 = 5.85 (6.0) 
by a Cabibbo-Marinari-Okawa algorithm with 8 hit pseudo heat bath algorithm for three 
577(2) subgroups. The acceptance rate is about 0.95 for both /3's. Each configuration is 
separated by 1000 sweeps after a thermalization of 6000 (22000) sweeps at (3 = 5.85 (6.0). 

The quark propagator G on a configuration given by 

D(K, n, m)G(m) = B(n) (4) 

m 

is constructed using a red/black minimal residual algorithm, taking periodic boundary con- 
ditions in all directions. We employ the point source at the origin B(n) = <5 n . 

The convergence criterion we take for the quark matrix inversion is that both of the 
following two conditions be satisfied: 

y/\R\ 2 /(3 xAxV) < 10~ 9 , (5) 

max niCiS {|-R CiS (n)/G CiS (n)|} < 0.03, (6) 

where \R\ is the norm of the residual vector R = B — D(K)G, V = L 3 x T is the lattice 
volume (L = 24 is the lattice size in the spatial directions and T = 54 is that in the temporal 
direction), and c and s are color and spin indices. The average number of iterations needed 
for the convergence is given in table 0. 

Selecting several configurations, we have solved exactly eq. [| within single precision to 
construct an exact hadron propagator and compared it with that obtained with the stopping 
conditions above. We find that the difference in a hadron propagator (for any particle at 
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any time slice) is at most one percent of the statistical error estimated using all (100 or 200) 
configurations. Therefore the error due to truncation of iterations is small enough and does 
not affect the following analyses and results. 

We use uTd for meson operators with T = 75 for ir, 27075 for n (fr), and 7, for p. For 
baryons, we use non-relativistic operators 

Nt = e abc j2K^ 'dfa 1 = 1,2 (7) 
A, = e abc £ Sf u1u)u c k I = ±3/2, ±1/2, (8) 

i,j,k 

where T3 is the third component of Pauli matrices and Si is the projection operator to 
J = 3/2, J z = I state. We also use anti-baryon operators obtained by replacing the upper 
components of the Dirac spinor in eqs. |7| and ^ with the lower components. 

We average zero momentum hadron propagators over all states with the same quantum 
numbers; three polarization states for the p meson and two (four) spin states for the nucleon 
(A). Then we average the propagators for particle and anti-particle: For mesons we average 
the propagator at t and that at T — t, for baryons we average the propagator for particle at 
t and that of anti-particle at T — t. In this work we only calculate the masses of hadrons 
composed of degenerate mass quarks. 

Statistical independence of hadron propagators on each configuration is investigated by 
the following two methods. 1) We divide the total propagators into bins of Nb successive ones 
and apply the single elimination jack-knife method to N con f/NB block-averaged propagators. 
We find that the errors in various quantities do not change significantly even if we change 
the bin size. Fig. [TJ shows typical results for the bin size dependence of the error in effective 
masses. 2) If configurations are independent, we expect that the error obtained for the set 
of iV configurations A(JV) behaves as 

A(N) ~ 1/Vn. (9) 

We check that this behavior is approximately satisfied using the propagators on the first iV 
configurations. Fig. ^ shows typical results for the N dependence of the error in effective 
masses. 



3 Hadron Masses 
3.1 Fitting procedure 

Ground state masses of hadrons are extracted by fitting hadron propagators G(t) to their 
asymptotic forms: 

G Q (t) = A{exp(-mt) ± exp(-m(T - £))} (10) 

for mesons and 

G (t) = Aexp(-mt) (11) 

for baryons. (We will discuss the masses of excited states later.) We perform least mean 
square fits taking account of time correlations minimizing x 2 defined by 

tmax 

X 2 = E {G(t)-G (t)}C- l (t,t'){G(t')-G (t')} (12) 

t f'—t 
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where C~ l (t,t') is the inverse of correlation matrix C(t,t') it min <t,t'< t max ). Errors are 
estimated by two methods. One is the single elimination jack-knife method taking account 
of the correlations among the propagators at different time separations. Another estimate 
of the error is obtained from the least mean square fit itself. Linear approximation to the 
fitting function around the minimum of x 2 gives a linear relation between the variance of fit 
parameters and the variance of propagator G(t) for the fitting range t = t min — t max . The 
relation leads to the error propagation rule which relates the correlation matrix C (t, t') to 
the error (and the correlation) on the fit parameters. We find that the errors obtained by 
the two methods are of the same order and that the error obtained by the jack-knife method 
is slightly (0% to at most 40%) larger than that by the least mean square fit. Hereafter we 
quote the former error for the sake of safety, unless otherwise stated. 



3.2 Fitting ranges and systematic error analyses 

In order to obtain a ground state mass, we have to choose carefully the fitting range t min 
- t max in such a way that the contamination from excited states is negligibly small. We 
fix tmax — T/2 in order to take into account the data at as large distances as possible. For 
the purpose of fixing t min , we make fits to a range t — T/2 varying t which is a candidate 
for t min . Then we investigate the t dependence of the fitted mass mf it and x 2 /df, df being 
the number of degrees of freedom, together with the t dependence of the effective mass m e ff 
defined by 

G(t)/G(t + 1) = G (t, m eff (t))/G {t + 1, m eff (t)). (13) 

We plot in Figs. ||, £| and ||, as examples, the results for x 2 /df, mf it and m e jj at f3 = 6.0, 
K = 0.155 for the pion, the p meson and the nucleon, respectively. Common features of the 
time slice dependences of x 2 /df, mf it and m e ff for all cases including the other cases which 
are not shown here can be summarized as follows. (Discussion on each particle together with 
a complete set of figures for effective masses will be given below.) 

1) When we increase to starting from a small value such as to = 4, x 2 /df decreases rapidly 
from a large value down to a value around 2.0 — 0.5 and stabilizes. We denote to where the 
stabilization starts as t x 2. The stabilized value of x 2 /df depends on particle, (3 and K. In 
table we give t x 2 and x 2 /df at t x 2. We note that t x 2 for lighter quarks are smaller than 
those for heavier quarks. From a point of view of the least mean square fit, t x 2 as well as 
any value of t > t x 2 are candidates for t min . 

2) Although m e ff(t) and m/^(t) almost stabilize around t ~ t x 2, a clear long plateau in 
m e ff is rarely seen and the data of m e jf frequently show large and slowly varying fluctuations 
at large time separations, as shown in the figures. If the fitting range is fixed case by case 
based on a short plateau of m e jj , this may lead to a sizable underestimate of statistical errors. 

3) The value of m e jj in many cases is still decreasing at t ~ t x 2. Similar phenomena 
are reported by the UKQCD collaboration [JEJ. Although probably the large statistical 



fluctuations mentioned above is a partial cause of this phenomenon, the possibility that 
excited states still contribute at t ~ t x 2 can not be excluded. It is difficult to clearly 
separate out the effects of excited states from the statistical fluctuations. 

From these considerations, we do not simply take t x 2 as t min . In order to remove the 
contamination from exited states as much as possible, we proceed in the following way: We 
take t m i n common to all K's for the mesons and for the baryons, respectively, at each /3, 
in order to avoid a subjective choice case by case. Therefore, we require t min > t x 2 for all 
K's. We further require that t m i n always lies in a plateau when a clear plateau is seen in the 
effective mass plot. In cases where two plateaus are seen (e.g. see Figs. ^ — 15|), we require 
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that t min is larger than the beginning point of the first plateau. We also pay attention to 
the consistency between the choices of t m i n at two /3's in such a way that the ratio of the 
values of t m i n is approximately equal to that of the lattice spacings at two /3's. Thus we 
have chosen t min = 12 (15) for mesons and t min = 13 (16) for baryons at (3 — 5.85 (6.0), 
respectively. The ratio of t min at (3 = 5.85 to that at (3 — 6.0 is approximately equal to the 
ratio of the lattice spacings a((3 = 5.85)/a(/3 = 6.0) ~ 1.2. 

In addition to statistical errors, we estimate the systematic error coming from uncertain- 
ties in the choice of fitting range 0, [l2]| . Varying to from t x 2 up to t min + 4, we estimate 
the upper (lower) bound for the systematic error by the difference between the maximum 
(minimum) value and the central value obtained from the fit with to = tmin- We take to 
only up to t min + 4, because, when t is larger than this value, data in the fitting range 
become too noisy. (For the A baryon at (3 = 6.0, we vary to up to t m i n + 3 because a fit with 
= tmin + 4 = 20 does not converge.) 

In this way we estimate the errors in ground state masses due to statistical fluctuations as 
well as those due to the possibly remaining contamination from excited states which cannot 
be properly taken into account by the standard least mean square fit with a fixed fitting 
range. Note that the data are consistent with the implicit assumption that the ground state 
dominates for t > t m i n when we take into account these systematic errors. Consistency of 
this assumption is also checked by a two-mass fit discussed in Sec. |[ 

3.3 Pion masses 

We show m e ff at (3 = 5.85 and (3 = 6.0 in Fig. || The pion effective mass has structure 
with the scale of the standard deviation even for t>t x 2\ In some cases m e jj(t) exhibits two- 
plateau structure or slow monotonic decrease. However, the magnitude of the fluctuation for 
the pion is much smaller than the other cases. The resulting systematic error is comparable 
to the statistical uncertainty. The results of the fits are given in table |3|. 

3.4 Rho meson masses 

Fitting to the p meson propagator is more problematic than the pion propagator. Because 
of this, we will discuss it at some length and compare the results with previous works. 

The p meson effective mass at (3 = 5.85 shown in Fig. 0-a exhibits a plateau for t > t x 2 = 
12 for the smallest two K's, while it exhibits peculiar behavior at large t for the largest three 
i^'s: m e ff(t) for t = 17 — 20 is larger than that for t = 12 — 16 and it drops abruptly at 
t = 21. We regard this behavior as due to statistical fluctuations. We find that fits to a 
range t = 12 — t max are stable for t max = 14 — 27. Therefore we choose t max = T/2 even for 
these cases. The results of the fits are summarized in table |j. The systematic error upper 
bound is 1 — 2 times larger than the statistical error for the largest three K's. 

Fig. 0-b shows the effective mass at /3 — 6.0. Except for the smallest K, m e ff(t) is 
decreasing at t ~ t x 2. Rate of the decrease becomes slow at t ~ 12 to exhibit a plateau 
for two or three time slices. The value of m e ff decreases further up to t ~ 17 to attain 
another plateau. The plateaus are not long enough to determine unambiguously the time 
slice where the contribution of excited states can be ignored. It should be emphasized again 
that x 2 /df are almost identical both for the fits with t m i n = 12 and t m i n = 17: 1.35 and 1.16 
for K = 0.1550, 1.20 and 1.13 for K = 0.1555 and 0.77 and 0.76 for K = 0.1563, respectively. 
See also Fig. [|. Therefore the value of x 2 does not give a guide to determine t min . The point 
tmin = 15 is located between the two pseudo-plateaus at t ~ 12 and t ~ 17. In table [| 
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are summarized the results for the fits with t min = 15 together with the systematic error. 
Reflecting the slow monotonic decrease of effective masses, the ratio of the systematic error 
to the statistical error is relatively large: the systematic error amounts to about twice the 
statistical error for the largest three i^'s. 

We notice a very intriguing fact that mf it by the correlated fits to a range from t = to to 
T/2 has a strong correlation with m e ff at t = t Q . A typical example is seen in Fig. |j. This 
holds for the other particles also. This means that the result of the fit to a range to — T/2 
is mainly determined from data at t ~ to- 
la our previous work H, we analyzed the same set of p meson propagators with un- 
corrected fits. Paying attention to the monotonic decrease of effective masses, we made 
two different fits to estimate the systematic error coming from uncertainties in the choice of 
fitting range. One is a fit to a range t ~ 9 — 11 at (3 — 5.85 (t ~ 12 — 15 at (3 = 6.0). We 
called the fit "pre-plateau fit". Another is a fit to a range t ~ 11 — t max at (3 = 5.85 (t ~ 
15 — t m ax at (3 — 6.0), which we called "plateau fit". The latter fitting ranges correspond 
approximately to those we adopt in this work. Because m e jj is decreasing, the p masses 
obtained from the correlated fits are systematically larger than those from the uncorrelated 
fits, due to the fact give in the preceding paragraph. The mass value obtained in this work 
is between that from the uncorrelated plateau fit and that from the pre-plateau fit. 

In table ^ we reproduce the results for the p meson masses at (3 = 6.0 for K = 0.155 
and 0.1563 together with those by the Ape collaboration || [7[j and the LANL group [II|. In 



1991, the Ape collaboration reported the result obtained on a 24 3 x 32 lattice with a multi- 
origin 7 3 cubic source j7|. Then we made simulations for the same spatial size with larger 
temporal extent [|J, 24 3 x 54, using the point source. For K = 0.155, the values of m e g at 
t ~ 10 are in close agreement with Ape's. Consequently the result 0.4280(33) obtained from 
the pre-plateau fit (t = 12 — 15) agreed with the Ape result 0.429(3) within one standard 
deviation. However, the result 0.4169(48) from the plateau fit (t = 15 — 27) was smaller 
by approximately twice the statistical error. We regarded the latter more reliable. At that 
time there was a report that the mass value appears to depend on the type of the source 



adopted |T3|. Therefore, in order to clarify whether the origin of the discrepancy between 
our result and the Ape result is due to the different type of source, we made calculation 
at K — 0.155 for 400 configurations ]5| using the point source, the wall source and the 
source adopted by the Ape collaboration. The results obtained from correlated fits for the 
three different sources have agreed with each other; 0.4201(29), 0.4228(19) and 0.4249(19) 
for the point source, the wall source and the multi-origin source, respectively. The recent 
result reported by the LANL group 0.422(3) JTTJ] is consistent with these numbers. It is 
probable that the slightly large value by the Ape collaboration is due to small temporal 
extent. The Ape collaboration has also made simulations using both the point source and 
the multi-cube source || with larger temporal size and smaller spatial size: 18 3 x 64. Their 
results 0.430(10) and 0.428(8) are consistent with other results within relatively large errors, 
although the central values are slightly higher than the results by other groups. The slightly 
large central values may be due to the small spatial size. For K = 0.1563, the results obtained 
from the correlated fit in this work is consistent with those by the Ape collaboration and 
the LANL group, albeit with large errors in the results. 

3.5 Baryon masses 

Fig. |] shows effective masses for the nucleon at (3 = 5.85 and (3 = 6.0. Decrease of m e jj at 
t ~ t x 2 is not conspicuous compared with the case of the p meson. However, we see two- 
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plateau structure for the cases of K — 0.1585 and 0.1595 at /3 — 5.85 and K = 0.155 (see also 
Fig. ||) and 0.1555 at (3 — 6.0. The choice t m i n = 13 (16) for (3 = 5.85 (6.0) corresponds to 
that we select the first (last) plateau as correct for the case where two plateaus are observed. 
Table [J] summarizes the results of the fits. 

For A, monotonic decrease of effective masses at t ~ t x 2 or two-plateau structure is seen 
for K = 0.1595 and 0.1605 at (3 = 5.85 and for K = 0.1550 and 0.1563 at /3 = 6.0. Effective 
mass plots are shown in Fig. |]. The results of the fits are summarized in table [7[ 

In table the baryon masses at (3 = 6.0 for K = 0.155 and 0.1563 together with those by 
the Ape collaboration and the LANL group are reproduced. The nucleon masses reported 
by the three groups agree within the statistical uncertainties. The A masses for K = 0.155 
are slightly scattered: Our result is higher than the LANL result by two standard deviations. 
However note that the values of the A mass obtained on 400 configurations (0.7054(95), 
0.7008(57) and 0.7128(191) for the point source, the wall source and the multi-origin source, 
respectively) are in good agreement with the LANL result. Therefore we think that the 
difference between the LANL result and our present result is due to statistical errors. 



3.6 Finite lattice effects 

The linear extension of the lattice in the spatial directions is 2.45 (2.03) fm at (3 = 5.85 
(6.0), when we use a" 1 = 1.93 (2.33) GeV determined from m p (see Sec. |5|). These values 
are much larger than twice the electromagnetic radius of the nucleon, 2x 0.82 fm. We also 
note that our results on the lattice with spatial volume 24 3 agree well with those on a lattice 



with 32 111, as discussed above. Therefore we do not take into account in this work finite 



lattice effects which are supposed to be small. 



3.7 Mass ratios 



The mass ratio m^/mp is plotted versus (m^/m-p) 2 in Fig. 10|. The values of the mass ratio 
are given in table [8]. The value of m^jm p at (3 = 6.0 is systematically smaller than that at 
(3 = 5.85, although the results at two /3's agree within the statistical uncertainty except for 
the case of the heaviest quark ((m n /m p ) 2 ~ 0.94). 



4 Excited State Masses 

In addition to the masses of ground states, we study the masses of first excited states for 
the p meson and the nucleon. To this end, we perform two-mass fits to the corresponding 
propagators varying t min . Our results for the p meson are shown in Fig. |TT] for (3 = 5.85, 



K = 0.1585, and in Fig. [12] for (3 = 6.0, K = 0.155. The results for the nucleon are given in 
Figs. [T3] and [T4] for (3 = 5.85 and 6.0, respectively. We find the following: 

1) x 2 /df is stable and small (~ 1 — 2) for t m i n > 4 (5) in the case of the p meson and 
for t min > 5 (6) in the case of the nucleon at f3 = 5.85 (6.0), respectively. 

2) When x 2 /df is small, the ground state masses ttiq from the two-mass fit are consistent 
with those from the one-mass fit within the errors, although the errors for mo from the 
two-mass fit become extremely large at large t m i n . 

3) Although x 2 /df is stable, the mass of the first excited state mi is in general quite 
unstable. For example, for the p meson at f3 = 5.85, the value of mi decreases from 1.5 for 
tmin = 3 to 0.6 for t min = 9 (cf. Fig. [□]). Similar behavior is also seen in the results for the 
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p meson at (3 — 6.0 (Fig. [12|) and the nucleon at f3 = 5.85 (Fig. |13|). The case of the nucleon 
at (3 = 6.0 is exceptional: mi is relatively stable (Fig. [TJ]). 

Under these circumstances, we select two t m j n 's which give m consistent with the result of 
the one-mass fit, under the condition that the errors are small. We then investigate whether 
the results for the excited state mass are consistent with the corresponding experimental 
values. 



In Figs. 15 and 16 are shown the first excited state masses of the p meson obtained 



from the fit with t min = 5 and 6 (8 and 9) versus 1/K at (3 = 5.85 (6.0), respectively. (A 
two-mass fit with t min = 9 for the largest K at (3 — 6.0 does not converge. Therefore the 
corresponding data is missing in the figure.) We give in the figures the experimental values 
for the masses of p(1450) and 0(1680) which are the first excited states of the vector mesons. 
The mass of 0(1680) is plotted at the third largest K, because this value of K corresponds 
to the strange quark mass as mentioned in Sec. |5.71 . Apparently the results for the excited 
state mass depend strongly on the value of t m i n . For quarks lighter than the strange quark, 
the excited state mass obtained with smaller t min is much larger than experiment, while 
that with larger one is consistent with experiment within large statistical errors. Therefore, 
although the value of m,\ is unstable, there exist two-mass fits to the p propagators which 
give both the ground state mass consistent with the one-mass fit and the first excited state 
mass consistent with experiment. 

Fig. [T7] shows the masses of excited state of the nucleon at (3 = 5.85 versus 1/K. The 
excited state masses obtained from the fit with t min = 7 are much smaller than those with 
tmin = 6. (A two-mass fit with t min = 7 for the largest K does not converge.) We expect 
that the mass difference between the ground state and the first excited state depends only 
weakly on the quark mass, because the mass difference for the spin 1/2 baryon satisfies this 
property. The mass difference for the nucleon is m7v(i44o) — ^tv(94o) = 500 MeV. The figure 
shows that the excited state masses with t min = 7 lie approximately 500 MeV higher than 
the ground state masses. Therefore there exist two-mass fits whose results do not contradict 
with experiment also for the nucleon at [3 — 5.85. 

In Fig. ^ we show the excited state masses of the nucleon at (3 = 6.0 with £ m j„ = 7. The 
masses of the first excited state lie much more than 500 MeV above the ground state masses. 
As mentioned before, two-mass fits for the nucleon at f3 — 6.0 are stable and therefore the 
values of the excited state mass do not change much even if we take other t min . When we 
recall that there exists a fit which gives a reasonable excited state mass at (3 — 5.85, this 
situation is puzzling. One possible origin for the heavy excited state mass at (3 = 6.0 is 
a finite size effect, because the physical volume is smaller at (3 = 6.0. There remains a 
possibility that when we simulate on a larger lattice, a two-mass fit with larger t min gives a 
value consistent with the nucleon excited state mass. 

There are several published data for the mass of excited states |5|, 0, In table 

|], we reproduce the results for the ratio of the excited state mass to the ground state mass 
selecting the quark mass corresponding approximately to the strange quark mass. For the p 
meson, except our results in this work with t min = 6 (9) at /3 — 5.85 (6.0) and the result for 
the wall source in ref. ||, the reported ratios are considerably larger than the corresponding 
experimental value m^ieso)/™^™) = 1-65. For the nucleon, the mass ratios reported by 
the Ape collaboration and the UKQCD collaboration are considerably larger than our result. 
One possible origin of the differences is due to the choice of fitting range. Because the two- 
mass fit is very unstable, we certainly have to employ a more efficient way to extract reliable 
values for the excited state masses. 



8 



5 Masses of Hadrons with Physical Light Quarks 



5.1 Extrapolation procedure 

Extrapolation of hadron masses to the chiral limit is done with the correlation being taken 
into account, among the masses at different values of hopping parameter. First we consider 
a least mean square fit to minimize 

X 2 = E {G(t,K)-Go(t,^)}C- 1 (t,K;t',ir'){G'(t , ,^ , )-G'o(t',^ , )} ) (14) 

t,t',K,K> 

where Go(t, K) = A(K)e~ m<yK ' >t is the fitting function to hadron propagator G(t, K) and C^ 1 
is the inverse of the full correlation matrix C(t, K;t'; K'). A linear approximation to the 
fitting function around the minimum of x 2 gives the relation between the error matrix E for 
fit parameters and correlation matrix C(t, K;t' , K') for propagators: 

S = (D T C- 1 D)- 1 , (15) 

where D is the Jacobian defined by 

D t , K ;A(K')MK') = [9G (t, K)/dA(K'), dG (t, K)/dm{K')\. (16) 



(D is diagonal with respect to K.) The full least mean square fit to minimize x m ec L- H 
is different from the set of least mean square fits for each K to minimize x 2 's in eq. 12: The 



masses and amplitudes obtained by the two methods are in general different. We take those 
obtained from the fits to each propagator for evaluation of the Jacobian.0 
For extrapolation, we minimize x 2 given by 

X 2 = J2{m(K) - f(K)} Z-\K, K') {m(K') - f(K')}, (17) 

K 

where the correlation matrix Ti(K, K') is the sub-matrix among the masses of the full error 
matrix S and f(K) is the fitting function. (For the pion, m(K) is replaced by m 2 (K) with 
appropriate replacement of S _1 (ii', K').) 

5.2 Linear extrapolation to the chiral limit 

We fit the data of the mass squared for the pion and the mass for the other hadrons at the 
largest three K's to a linear function of 1/K; f(K) = a® + a%/K. We find that quality of 
the linear fit is good in the sense that x 2 /df < 2 (df = 1 in this case) and therefore we do 
not study in this work the effects of possible chiral logarithms [f24| , p5| . We summarize the 
fit parameters together with x 2 /df in table |TU[ The linear extrapolations of hadron masses 
at (3 = 5.85 and 6.0 are shown in Figs. [19| and BUI respectively. 



In table |TT] we summarize the results for the critical hopping parameter K c and the masses 
at K c together with the errors estimated by the least mean square fit and those by the jack- 
knife method. We find that the error estimated by the jack-knife method is larger than that 
by the least mean square fit except for K c at (3 = 6.0. We take the error obtained by the 
jack-knife method as our estimate of the statistical uncertainty, unless otherwise stated. 



1 We have checked that the error matrix thus obtained is very close to that obtained using the Jacobian 
at the absolute minimum of eq. Consequently the difference in the extrapolated values obtained using 
two error matrices is at most 5% of their statistical uncertainties. 
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5.3 Systematic error analyses 

We first estimate the systematic error on the masses in the chiral limit coming from uncer- 
tainties in the choice of fitting range for extracting the ground state mass at each K. To 
this end, we repeat linear extrapolations of the masses obtained from the fits to a range 
t — T/2, varying t (common to all ICs) from max K {t x 2(K)} to t min ± 4. We find that 
quality of the linear fits depends on the choice of to: X 2 /df are considerably large for some 
choices of t Q . We adopt the condition x 2 /df < 2 for the linear fit to be accepted. We take 
the difference between the fitted mass value and the maximum/minimum mass value under 
the condition x 2 /df < 2 as our estimate of the systematic upper/lower error. We call the 
systematic error thus obtained the fit-range systematic error. 

Data at the fourth largest K slightly deviates from the linear fit. In order to estimate 
the systematic error which comes from the choice of fitting function, we make a quadratic 
fit {f(K) = do ± di/K ± a 2 /K 2 ) to the largest four K's, varying to in the range used for the 
estimate of the fit-range systematic error. We estimate the systematic error by the difference 
between the maximum/minimum value with x 2 /df < 2 and that of the linear fits. We call 
the systematic error thus obtained the fit-func systematic error. 

5.4 Pion mass extrapolation and K c 

Pion masses squared are fitted to a linear function of l/K to obtain the critical hopping 
parameter. The value of x 2 /df is 0.56 (1.1) for the fit (t min = 12 (15)) at (3 = 5.85 (6.0). 
The fit-range systematic errors are estimated from the fits with to — 8 — 16 at (3 — 5.85 
and 10 — 19 at (3 = 6.0. All the fits give x 2 /df < 2. The upper (lower) bound comes from 
the fit with t = 11 (14) with X 2 /df of 0.36 (0.04) for (3 = 5.85 and from the fit with t = 12 
(19) with x /df of 0.44 (0.96) for (3 = 6.0. 

For data at (3 = 5.85, no quadratic fits with t = 8 — 16 give x 2 /df < 2. On the other 
hand, quadratic fits to data at (3 = 6.0 with t = 13 — 19 give x 2 /df < 2. Because m 2 , is 
a concave function of l/K when the data at the fourth largest K is included, K c obtained 
from the quadratic fit is larger than that from the linear fit. 

The values of K c 's together with the fit-range systematic error and the fit-func systematic 
error are given by 

stat. sys. (fit-range) sys. (fit-func.) 

(3 = 5.85 K c = 0.161624 ±0.000033 +0.000001 -0.000025 
P = 6.00 K c = 0.157096 ±0.000028 ±0.000033 -0.000009 ±0.000109 

The fit-range systematic error is comparable to the statistical uncertainty. 

The result for K c at (3 — 6.0 agrees well with that in ref. 0. Although it is slightly 
smaller than the LANL result 0.15714(1) [|ll|], we conclude that our result is consistent with 
theirs within the sum of the statistical error and the fit-range systematic error. 

In this work, we do not distinguish the physical point where m^jm p takes its experimental 
value from the critical point where the pion mass vanishes, because we find that physical 
quantities at the two points differ only at most 30% of their statistical errors. 

5.5 Rho meson mass extrapolation and lattice spacing 

A linear fit to the p meson masses (with t m i n = 12 (15)) at the largest three i^'s gives x 2 /df 
of 1.8 (1.2) for (3 = 5.85 (6.0). Therefore the linear fit is acceptable. 
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However, we find that quality of the linear fit strongly depends on the choice of fitting 
range. See Figs. |2"T| and |2"2]. In table we summarize x 2 /df ', m p (K c ) and the inverse lattice 
spacing defined by a -1 = 0.77GeV /m p (K c ) versus £q. 

We also make a quadratic fit to the data at the largest four i^'s to estimate the systematic 



error due to the choice of fitting function. Table summarizes the results of the quadratic 
fits versus t . 

The method to estimate the systematic error is the same as that adopted for the pion. 
Our final results for m p are 

stat. sys. (fit-range) sys.(fit-func) 

P = 5.85 m p {K c ) = 0.400 ±0.021 +0.008 -0.027 +0.0 -0.013 
P = 6.00 m p {Kc) = 0.331 ±0.011 +0.018 -0.020 +0.0 -0.008 

The value of m p (K c ) at (3 — 6.0 agrees well with the Ape result 0.3332(75) and the LANL 
result 0.3328(106). The values of m p (iT c )'s are translated to the lattice spacing as 

stat. sys. (fit-range) sys.(fit-func) 

= 5.85 a- 1 = 1.93 ±0.10 +0.14 -0.04 +0.08 -0.0 GeV 

/5 = 6.00 a' 1 = 2.33 ±0.08 +0.15 -0.12 +0.06 -0.0 GeV 

Although the statistical error on a -1 is several percent, we notice that the systematic error 
is much larger. Summing up both the statistical and systematic errors, we find that a -1 can 
be as large as 2.25 GeV (2.62 GeV) at P = 5.85 (6.0) and as small as 1.79 GeV (2.13 GeV). 

In analyses of the systematic errors above, we have taken t common to all K's. However, 
it is not necessary to restrict ourselves to take a common value of to, because the time slice at 
which the contribution of excited states becomes negligible can depend on the quark mass. 
We make linear fits to all possible combinations of the p masses at the largest three K's, 
varying t separately for each K from t x 2 to 18. Fig. ^ shows a -1 at {3 = 6.0 versus x 2 /df- 
We see that there are linear fits with small x 2 /df which give both large and small a -1 . The 
value of a -1 scatters approximately from 2.15 GeV to 2.65 GeV. This upper value as well as 
the lower value are consistent with those obtained above with the systematic errors included. 

We estimate the value of J defined by my |23j from the linear fits discussed above: 



stat. sys. (fit-range) 

P = 5.85 J = 0.420 ±0.049 +0.028 -0.024 
p = 6.00 J = 0.395 ±0.026 +0.026 -0.026 



The value of J at P — 6.0 is smaller than the experimental value 0.48(2) even when we 
include the systematic errors. 

5.6 Nucleon and A masses 

Both linear fits and quadratic fits are made to the masses of the nucleon and the A baryon 
by the same method as for the p meson. Results of the linear fits versus the fit-range are 



summarized in tables and The fit with t m i n = 13 (16) at P — 5.85 (6.0), which is 
adopted in this work, gives a small x 2 /df = 0.37 (0.05). For the nucleon, quality of the 
linear fits is good for almost all values of to i n the sense that x 2 /df are approximately less 
than 2, except for the fit with to — 9 at P — 5.85. This feature is different from that for the 
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p meson. Quality of the fits to the A masses at /3 — 5.85 is good for t < 13 including our 
choice t m i n = 13 and that at (3 — 6.0 is good for all to except for t = 13. 
Results with various errors are give by 

stat. sys. (fit-range) sys(fit-func) 

(3 = 5.85 m N {K c ) = 0.589 ±0.036 +0.018 -0.058 +0.0 -0.018 
(3 = 6.00 mjviifj = 0.462 ±0.024 +0.020 -0.009 +0.0 -0.007 



stat. sys. (fit-range) sys.(fit-func) 

(3 = 5.85 m A (K c ) = 0.664 ±0.063 +0.034 -0.0 +0.0 -0.031 
f3 = 6.00 m A (K c ) = 0.605 ±0.033 +0.041 -0.011 +0.016 -0.007 



The value of the nucleon mass in the chiral limit at (3 = 6.0 lies between the LANL result 
0.482(13) and the Ape result 0.432(15). For the A masses, results by the three groups agree 
well with each other, albeit with large errors; the LANL result is 0.590(30) and the Ape 
result 0.58(3). The LANL results are those at the physical point where m n /m p takes its 
experimental value. 

These results are translated to the masses in physical units using the value of a -1 obtained 
from m p . The systematic error on the lattice spacing is not taken into account for the 
estimate of the systematic error on the baryon masses. Results read 











stat. 


sys.(fit- 


range) 


sys.(fit-func) 




p = 


5.85 


m N = 


1.135 


±0.088 


+0.034 


-0.112 


+0.0 


-0.034 


GeV 


p = 


6.00 


m N = 


1.076 


±0.060 


+0.047 


-0.020 


+0.0 


-0.017 


GeV 


p = 


5.85 


m A = 


1.279 


±0.136 
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6.00 
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+0.096 


-0.026 


+0.038 


-0.015 


GeV 



The central value of the nucleon mass at (3 = 6.0 (5.85) is larger than its experimental value 
by about 15% (20%) and that of the A mass by about 15% (4%): The errors amount to 
twice the statistical errors except for the A baryon at (3 = 5.85. The systematic errors 
are comparable with the statistical errors (3 - - 13%). Even when the systematic errors 
are included, the baryon masses at (3 — 6.0 do not agree with experiment. Our data are 
consistent with the GF11 data |I| at finite lattice spacing, within statistical errors. In order 
to take the continuum limit of our results, we need data for a wider range of (3 with statistical 
and systematic errors much reduced. 



5.7 Masses of strange hadrons 

The hopping parameters for the strange quark which are estimated from the experimental 
value of rax/ m p turn out to be K s = 0.1588 and 0.1550 at (3 = 5.85 and 6.0, respectively. 
Note that they are identical or almost identical to the third largest hopping parameter 
K = 0.1585 and 0.1550 which we have chosen in such a way that they approximately 
correspond to the strange quark. The masses of Q~ estimated at K = K s are 1.696(92) GeV 
and 1.693(57) GeV at f3 = 5.85 and 6.0, respectively (statistical errors only). They are in 
good agreement with the experimental value 1.672 GeV. The masses of the vector meson at 
K = K s are 998(45) MeV and 986(26) MeV at p = 5.85 and 6.0, respectively, which equal 
the 4> meson mass 1019 MeV within about one standard deviation. As is well known, there 
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are ambiguities in determination of the hopping parameter for the strange quark. When the 
hopping parameters for the strange quark mass are alternatively determined from m^/m p , 
they are equal to 0.1585 and 0.1547. The results for the Q~ mass at these hopping parameters 
are consistent with those above within one standard deviation. 



6 Meson Decay Constants 

6.1 Vector meson decay constants 

We evaluate vector meson decay constants defined by 

(0\(u lt d) cont \V(p = 0)) = e t F v m v , (18) 

where e, and my are the polarization vector and the mass of the vector meson, respectively, 
and (u'jid) cont ' is the vector current in the continuum limit. The experimental value for the 
p meson is F p = 216(5) MeV. (This Fy is related to fy 1 by fy 1 = Fy/my.) 

The expectation value of the local lattice current {u^id) latt " between the vacuum and the 
vector meson is related to the continuum one by the relation 

(0\(u ll d) cont \V(p = 0)) = Z K Z v (0\(uy i d) latt -\V(p = 0)). (19) 

The coefficient Zk is a scale factor for the difference between the continuum and lattice 
normalizations of the quark field. The renormalization constant Zy is the ratio of the 
conserved lattice current to the local current, which can be estimated by perturbation theory 
or numerical simulations. We test the following three possible choices of Zk and Zy: 

1. those in naive perturbation theory: Zk = 2K and Zy = 1 — 0.174 g 2 [|T6f . 



2. those in tadpole improved perturbation theory: Zk = (1 — 3K/AK C ) |L7| and Zy = 
1 — 0.82a^r(l/a) |18| ( a Jds( 7! '/ a ) = 9jYs( 1T / a )/^ 71 * s determined by the relation 
l/ g 2—(n/a) = Ti(U P /3)/g 2 + 0.02461 pf ]l|. We then determine a m (l/a) using 
the two loop renormalization group equation.), 



3. Monte Carlo estimate of Z v = 0.51 [| (0.57 ||) at = 5.85 (6.0) with Z K = 2K. 
(Data for Zy at (3 = 5.85 are given in table [16]. Because the results for Zy are 
independent of the quark mass in the range we investigate, we use the averaged value.) 
The error on Zy is ignored in the following. 

We abbreviate the decay constants obtained using the above three renormalization constants 
as Fy T Fy P and Fy IC , respectively. 

The statistical error is obtained by the jack-knife method. The systematic error is esti- 
mated varying to as in the case of mass calculation. The range of to is the same as that for 



the p mass. In table [17] we summarize the results for the decay constants at each K. We 
quote the error only for Fy P , because the errors for the others can be easily obtained from 
that for Fy P by multiplying the ratio of Z-factors. 

Fig. |24| shows Fy/my versus (mps/my) 2 together with the corresponding experimental 



values for p, u, <p and J/ip. Note that we can compare the numerical results with the 
experimental values for <ft and J/ip without extrapolation. The values with Fy IC at two /3's 
remarkably agree with each other. Furthermore they agree well with the experimental values 
for cf) and J/ip. This implies that scaling violation in Fy IC is small. On the other hand, we 
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find sizable scaling violation in F PT and F pp . They are off the experimental values for <p and 
J/ip by 40 — 100%. We find that Ftf c /m v , s at (3 = 6.0 agree well with the Ape data 0, g. 



In fig. |25] we depict the values of F$jm p versus m p a together with the GF11 result pT 
The values of the hopping parameter for the strange quark are given in Sec. |5.7| . Note that 
the values of F^ IC /m p agree with experiment already at m p a = 0.33 — 0.40 within 1 — 2 
standard deviations. The values of Fj p /m p are consistent with the GF11 result, although 
the central values are about la higher than the GF11 data. They are off the experimental 
value by 30 — 40% at these values of m p a. Linear extrapolation of our data to zero lattice 
spacing is consistent with experiment. 

The value of Fy in the chiral limit is obtained from a linear fit in terms of 1/ K in a similar 
way to that made for hadron mass extrapolation. We first calculate the correlation matrix 
K') for Fy{K) from the error matrix £ for the mass and amplitude (eq. |15|) using the 
error propagation rule and then minimize \ 2 ■ A linear fit to the data at the largest three 
K'a gives a reasonable X 2 /df: X 2 /df = 0.04 (0.38) for F$ p , 0.09 (0.44) for F PT and F^ IC at 
(3 = 5.85 (6.0), respectively. Fig. ^6] shows F v as functions of the quark mass together with 
the fitting functions. 

The method to estimate the systematic error due to the choice of fitting range is similar 
to that for hadron masses at K c . The results of the linear fit for various fitting ranges are 
given in table 18. Our final results for F p read 
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0.111 


±0.008 


±0.016 
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F TP = 
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±10 


±37 


-40 


MeV 


pMC _ 


0.0944 


±0.0064 


±0.010 


-0.014 




pMC _ 


220 


±8 


±24 


-33 


MeV 



The values of F PT can be obtained from Fp IC by multiplying Z PT /Zy IC = 1.61 (1.45) at 
(3 = 5.85 (6.0). We show the values of F p /m p in Fig. |27|. It should be noted that the values 
of Fy C in the chiral limit at two /3's are consistent with the experimental value of F p . We 
find that our values of Fj p jm p are consistent with the GF11 result ||21|| , albeit the central 



values being roughly la lower than the GF11 data; this tendency is opposite to the case 
of the meson. We note that linear extrapolation of our data for F p TP /m p to zero lattice 
spacing is again consistent with experiment. 

6.2 Pseudo scalar meson decay constants 

The pseudo scalar meson decay constant is defined by 

(0\(u l0l5 dy ont -\P(p = 0)) = V2m PS fp S . (20) 

The experimental value is f n = 93 MeV. We investigate three cases of renormalization 
constants as in the case of Fy\ 1) = 1 — 0.133 g 2 in naive perturbation theory [l(| with 
Z K = 2K, 2) Z A = l-0.31o^(l/a) |B] with Z K = (1-3K/AK C ) [TJ in tadpole improved 



perturbation theory, and 3) Z& = 0.69 p0| at (3 = 6.0 as a nonperturbative evaluation with 
Zk = 2K. (Corresponding Z A at (3 = 5.85 is not known.) 
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We derive fps from a fit to the it propagator. The value of t m i n is chosen to be the same 
as that for tt. The pion mass from the tt propagator is given in table O. Although the mass 



obtained is 1 — 2 standard deviations smaller than that from the n propagator, they are 
consistent with each other if we take account of the systematic error. The decay constant 
at each K is given in table 20|. Our data for /pj at (3 = 6.0 and K = 0.155,0.1563 are 



consistent with the Ape results |7|, ||. Fig. |28| shows fps/ m v versus {mps/my) 2 together 
with the corresponding experimental values for tt and K and the upper bound for the D 
meson. Contrary to the case of the vector meson, fpZf differs from the experimental value 
for the K meson by a factor of about 1.2. There is a possibility that the lattice size 10 3 x 20 
is not large enough to suppress finite lattice size effects in the Monte Carlo evaluation of Za- 
We think we have to calculate nonperturbatively Za both at (3 — 5.85 and 6.0 on a larger 
lattice in order to clarify the reason of the discrepancy. 

In fig. |25| we show the values of f^/rrip versus m p a together with the GF11 result |pT| . 
The values of fx are evaluated at the hopping parameter given by 2/(l/K c ± 1/K S ). The 
values of fjf '/m p are consistent with the GF11 result, albeit with larger errors in our results. 
Our data at finite lattice spacing are also consistent with experiment. 

The extrapolation to the chiral limit is problematic. We find that neither of the linear 
fit to the data at the three largest K's nor the quadratic fit to the data at the four largest 
K's gives x 2 /df small enough: For fpg at (3 = 5.85 (6.0), x 2 /df = 9.1 (6.7) for the linear 
fit and 9.9 (7.4) for the quadratic fit, respectively. Fits to fp$ are similar. In Fig. ^ are 
shown fps versus the quark mass together with the linear fits. The data at the largest K is 
much below the fitting lines. Even if we change t min , x 2 Jdf does not reduce much. In Fig. 
|30| we show y 2 jdf together with the result for fpjj at (3 — 6.0 versus t min . Although x 2 /df is 
large, the results of the fits are very stable. Therefore we quote the decay constant obtained 
by the linear extrapolation of the data with t min = 12 (15) at (3 = 5.85 (6.0) as the central 
value of the decay constant. We estimate the systematic errors similarly as in the previous 
cases with t = t x 2 — 14 (16) for (3 = 5.85 (6.0). 

These analyses give 
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The values of f n obtained with the tadpole improved renormalization constants are con- 
sistent with the experimental value within the statistical errors (see Fig. 1271). That with 
the MC renormalization constant is also consistent with experiment if we take account of 
the (small) systematic error. However, we should take these numbers with caution, because 
X 2 /df for the extrapolation is not small enough as mentioned above. Note that the decay 
constants in the chiral limit are consistent with the GF11 data pl |, although the errors in 
our results are considerably larger. 
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7 Conclusions and Discussion 



In analyses of numerical simulations toward high precision determination of light hadron 
masses, one first encounters the problem of fitting range for hadron propagators. We find 
that effective masses of hadrons in general do not exhibit clear plateaus, although statis- 
tics is relatively high (the number of configurations is 100 (200) at /3 — 5.85 (6.0)). The 
correlated x 2 fits do not determine unambiguously the time slice beyond which the ground 
state dominates. We also notice a very intriguing fact that mfu by the correlated fits to a 
range from t — t has a strong correlation with m e g at t — t . Varying systematically the 
fitting range, we estimate systematic errors in hadron masses due to statistical fluctuations 
as well as due to the contamination from excited states, which cannot be properly taken into 
account by the standard least mean square fit with a fixed fitting range. We find that the 
systematic errors for the hadron masses with quarks lighter than the strange quark amount 
to 1 — 2 times the statistical errors. 

When the lattice scale is fixed from the p meson mass, the masses of the Q~ baryon and 
the meson at two /3's agree with experiment within about one standard deviation. On 
the other hand, the central value of the nucleon mass at /3 — 6.0 (5.85) is larger than its 
experimental value by about 15% (20%) and that of the A mass by about 15% (4%): Even 
when the systematic errors are included, the baryon masses at (5 — 6.0 do not agree with 
experiment. In order to take the continuum limit of the nucleon mass and the A mass, we 
need data for a wider range of (3 with statistical and systematic errors much reduced. For the 
masses of excited states of the p meson and the nucleon, there exist two-mass fits which do 
not contradict with experiment, except for the case of the nucleon at f3 — 6.0. Although this 
does not necessarily imply that the excited state masses appear consistent with experiment 
because two-mass fits are very unstable, the existence of such a fit consistent with experiment 
encourages us to perform more works in this direction. 

Determination of meson decay constants is usually accompanied by uncertainties of renor- 
malization constants. One can in principle employ any renormalization constant such as that 
determined by naive perturbation theory or tadpole improved perturbation theory. We have 
indeed shown that when we use renormalization constants given by tadpole improved per- 
turbation theory, although the decay constants for the 0, p, K and 7r mesons are in general 
off experiment at finite lattice spacing, for example, by 30 — 40% at m p a = 0.33 — 0.40 in 
the case of the F^, they approach in the continuum limit toward values consistent with the 
experimental values. 

It is, however, desirable to employ a renormalization constant which gives weak a de- 
pendence for the decay constants. We have shown that when we use the renormalization 
constants determined by Monte Carlo simulations, the vector meson decay constants at 
two /3's remarkably agree with each other and reproduce the experimental values within 
the errors for a wide range of the quark mass with the chiral limit included. This implies 
a strong advantage to apply renormalization constants determined nonperturbatively. For 
pseudo-sealer mesons, however, we find that although the decay constant fpg" in the chiral 
limit agrees with the experimental value of f n albeit with large errors, it differs from the 
experimental value of fx by about 20% at m p a = 0.33. This discrepancy might be due to 
systematic errors in the numerical calculation of Za- These results imply the importance of 
more systematic nonperturbative determination of the renormalization constants for various 
meson decays. 

Numerical simulations are performed under the QCDPAX project which is supported 
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by the Grants-in-Aid of Ministry of Education, Science and Culture (Nos. 62060001 and 
02402003). Analyses of data are also supported in part by the Grants-in-Aid of Ministry of 
Education, Science and Culture (Nos. 07NP0401, 07640375 and 07640376). 

Note Added 

After this work was completed, three groups have reported results of high statistics studies 
of the hadron spectrum []2B], ^7], ^] at /3 = 6.0. Their results are consistent with ours. 
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(3 = 5.85 


(3 = 6.0 


m n /m 


K 


-f-f- 1 "i" PT CI "f" 1 r\v\ 


K 


-Li- 1 \ pr cif inTi 
~f~f-L l)CI CL L ivjll 


0.1440 


80 ± 3 


0.1450 


90 ± 3 


97 


0.1540 


160 ± 10 


0.1520 


160 ± 10 


0.87 


0.1585 


420 ± 45 


0.1550 


380 ± 40 


0.70 


0.1595 


610 ± 75 


0.1555 


430 ± 45 


0.64 


0.1605 


1850 ± 410 


0.1563 


1110 ± 170 


0.52 



Table 1: Hopping parameters and average number of iterations used to solve quark propa- 
gators. Approximate values for m n /m p are also given. Table || contains precise values for 



(3 = 5.85 


K 


7T 


P 


N 


A 






X 2 /df 


f x 2 


X 2 /df 


l x 2 


X 2 /df 


f x 2 


X 2 /df 


0.1440 


12 


0.98 


12 


1.32 


11 


1.33 


11 


1.49 


0.1540 


10 


0.90 


12 


1.01 


11 


1.61 


11 


1.61 


0.1585 


8 


0.72 


8 


2.04 


9 


1.36 


11 


1.18 


0.1595 


8 


0.45 


8 


1.73 


9 


1.13 


11 


1.07 


0.1605 


8 


0.46 


8 


1.20 


7 


1.56 


9 


1.30 


= 6.0 


K 


71 


P 


N 


A 




t x 2 


X 2 /df 


t x 2 


xVdf 


t x 2 


XVdf 


t x 2 


xVdf 


0.1450 


15 


0.55 


15 


1.02 


15 


0.36 


15 


0.57 


0.1520 


12 


1.26 


13 


0.71 


15 


0.38 


15 


0.56 


0.1550 


10 


1.39 


11 


1.42 


12 


0.41 


12 


0.94 


0.1555 


10 


1.35 


10 


1.32 


12 


0.64 


12 


1.22 


0.1563 


9 


1.54 


9 


0.95 


10 


1.21 


11 


1.11 



Table 2: t x 2 and x 2 /df at t x 2. See the text for details. 
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(3 = 5.85 


(3 = 6.0 


K m n x 2 /df 


K m w x 2 /df 


0.1440 1.0293(12)^ 13.7/14 
0.1540 0.6122(11)11 10.9/14 
0.1585 0.3761(12)1| 7.4/14 
0.1595 0.3088(14)3 5.8/14 
0.1605 0.2226(21)3° 6.0/14 


0.1450 0.8069(7)1^ 6.1/11 
0.1520 0.4772(9)3 11.5/11 
0.1550 0.2967(15)3 8 14.1/11 
0.1555 0.2588(16)3 8 17.0/11 
0.1563 0.1847(27)lf 20.9/11 



Table 3: Pion masses in lattice units. In parentheses are errors estimated by the jack-knife 
method. Errors given in the form 3o««er are f° r the fitting range dependent upper/lower 
bound. 



(3 = 5.85 


(3 = 6.0 


K m p x 2 /df 


K m p x 2 /df 


0.1440 1.0598(15)3 18.5/14 
0.1540 0.6931(27)1?! 14.2/14 
0.1585 0.5294(69)1^ 28.9/14 
0.1595 0.4856(96)1}^ 23.2/14 
0.1605 0.434(20)1^ 14.6/14 


0.1450 0.8370(9)3 11.2/11 
0.1520 0.5486(15)1?! 6.6/11 
0.1550 0.4218(42)3| 14.3/11 
0.1555 0.3982(61)3f 12.4/11 
0.1563 0.353(15)1?? 7.6/11 



Table 4: The same as table [3] for the p meson. 



AT=0.155 



This work 24 3 x 54 
Ape 24 3 x 32 [§ 
Ape 18 3 x 64 || smear 
local 

LANL 32 3 x 64 |1| 



7T p 

0.2967(15) 0.4218(42) 

0.298(2) 0.429(3) 

0.297(2) 0.430(10) 

0.297(2) 0.428(8) 

0.297(1) 0.422(3) 



N 

0.6440(85) 
0.647(6) 



A 

0.728(11) 
0.745(15) 



0.641(4) 0.706(8) 



7^=0.1563 



7T 



This work 24 3 x 54 
Ape 24 3 x 32 [§ 
LANL 32 3 x 64 " 



P 



11 



0.1847(27) 0.353(15) 
0.184(3) 0.377(8) 
0.185(1) 0.363(9) 



N 

0.536(30) 
0.522(14) 
0.540(12) 



A 

0.670(53) 
0.636(45) 
0.631(27) 



Table 5: Comparison of hadron masses in lattice units at /3 — 6.0, K = 0.155 and 0.1563. 
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fj — o.oo 


r — r n 


K m N x 2 /df 


K m N x 2 /df 


0.1440 1.6961(50)^ 8 15.8/13 
0.1540 1.1060(55)1^ 22.3/13 
0.1585 0.815(13)1^ 17.5/13 
0.1595 0.744(17)1^ 18.1/13 
0.1605 0.683(48)1^ 23.4/13 


0.1450 1.3225(28)1^ 3.8/10 
0.1520 0.8669(49)1^ 9 4.2/10 
0.1550 0.6440(85)!?! 3.8/10 
0.1555 0.6007(109)lf 6.2/10 
0.1563 0.536(30)1^ 15.7/10 



Table 6: Nucleoli masses in lattice units. In parentheses are errors estimated by the jack-knife 
method. Errors given in the form 1"^^ are for the fitting range dependent upper/lower 
bound. 



P = 5.85 


P = 6.0 


K m A x 2 /df 


K m A x 2 /df 


0.1440 1.7124(57)1^ 17.7/13 
0.1540 1.1629(67)1^ 20.3/13 
0.1585 0.9011(153)!^ 16.5/13 
0.1595 0.825(21)1^ 15.1/13 
0.1605 0.755(53)1^ 19.4/13 


0.1450 1.3404(29)!f 6.2/10 
0.1520 0.9112(41)1^ 6.0/10 
0.1550 0.7278(109)1J 88 12.1/10 
0.1555 0.7001(159)1?^ 15.3/10 
0.1563 0.670(53)1^ 9.0/10 



Table 7: The same as table |6| for the A baryon. 



P = 5.85 


P = 6.0 


K m^/rrip m^jm p 


K m^/rrip m^jrap 


0.1440 0.9712(8) 1.6004(45) 
0.1540 0.8833(32) 1.5956(82) 
0.1585 0.7104(90) 1.540(29) 
0.1595 0.636(12) 1.531(42) 
0.1605 0.513(25) 1.57(12) 


0.1450 0.9641(5) 1.5801(25) 
0.1520 0.8699(21) 1.5802(79) 
0.1550 0.7033(69) 1.527(21) 
0.1555 0.650(10) 1.509(31) 
0.1563 0.523(23) 1.52(10) 



Table 8: Mass ratios m n /m p and m N /m p . The errors quoted are statistical only and are 
estimated by the jack-knife method. 
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p meson 


nucleon 




comment 


ratio 


comment 


ratio 


This work 


5.85 


train 5 


2.47(16) 


t ■ — 6 

b ram u 


1.64(12) 






train 6 


1.87(24) 


t ■ — 7 

b ram 1 


1.29(10) 




6.0 


train 8 


2.21(27) 


t ■ — 7 

b ram 1 


1.81(10) 






train 9 


1.58(26) 






Ape 


6.0 




2.13(21) 




2.13(4) 


UKQCD H 


6.2 




2.53(16) 




2.01(16) 


APE§ 


6.3 




1.93(10) 




1.93(12) 


UKQCD @ 


6.2 Clover 




2.23(14) 






QCDPAX § 


6.0 


point 


1.99(15) 


point 


1.55(20) 






wall 


1.70(26) 


wall 


1.47(21) 


experimental value 






1.65 







Table 9: Ratios of the excited state mass to the ground state mass. We have taken the quark 
mass corresponding approximately to the strange quark mass. 





(3 = 5.85 


P = 6.0 




a a x x 2 /df 


a ax x 2 /df 


m p 

m N 

m A 


-7.18(4) 1.16(1) 0.56 
-6.16(37) 1.06(6) 1.76 
-10.87(51) 1.85(8) 0.37 
-11.37(80) 1.95(13) 0.05 


-6.51(6) 1.02(1) 1.06 
-6.50(39) 1.07(6) 1.20 
-12.97(79) 2.11(12) 0.05 
-8.4(1.3) 1.42(21) 0.29 



Table 10: Fit parameters of the linear fits to the masses at the largest three iT's. Errors on 
oq and a± are those from least mean square fits. 



P = 5.85 




value 


err-lms 


err- jack 


jack/lms 


K c 


0.161624 


0.000027 


0.000033 


1.2 


m p (K c ) 


0.400 


0.010 


0.021 


2.1 


m N (K c ) 


0.589 


0.014 


0.036 


2.6 


m A (K c ) 


0.664 


0.022 


0.063 


2.9 


P = 6.0 




value 


err-lms 


err- jack 


jack/lms 


K c 


0.157096 


0.000038 


0.000028 


0.7 


m p (K c ) 


0.3309 


0.0080 


0.0114 


1.4 


m N (K c ) 


0.462 


0.015 


0.024 


1.6 


m A (K c ) 


0.605 


0.025 


0.033 


1.3 



Table 11: Values of K c and masses extrapolated to K c determined from the linear fits to the 
data at the largest three K's. Errors obtained by least mean square fits (err-lms) and those 
by the jack-knife method (err-jack) together with their ratios (jack/lms) are also given. 
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P = 


5.85 






P 


= 6.0 




to 


m p 


a- 1 


X 2 /df 


to 


m p 


a- 1 


X 2 /df 


8 


0.4359 


1.766 


36.90 


11 


0.3525 


2.184 


2.51 


9 


0.4213 


1.828 


11.17 


12 


0.3476 


2.215 


2.54 


10 


0.4196 


1.835 


14.80 


13 


0.3425 


2.248 


0.77 


11 


0.4045 


1.904 


9.30 


14 


0.3409 


2.259 


3.94 


12 


0.3998 


1.926 


1.76 


15 


0.3309 


2.327 


1.20 


13 


0.3892 


1.978 


4.39 


16 


0.3248 


2.370 


0.17 


14 


0.4081 


1.887 


1.36 


17 


0.3188 


2.416 


0.58 


15 


0.3794 


2.030 


0.43 


18 


0.3112 


2.474 


0.03 


16 


0.3728 


2.066 


0.11 


19 


0.3206 


2.401 


0.00 



Table 12: Results of the linear fits to the p meson masses versus t . The inverse lattice 
spacing is defined by a" 1 = 0. 77GeV /m p (K c ). 





P = 


5.85 






P 


= 6.0 




to 


m p 


a- 1 


X 2 /df 


to 


m p 


a- 1 


X 2 /df 


12 


0.3881 


1.984 


1.22 


13 


0.3413 


2.256 


0.79 


13 


0.3767 


2.044 


4.02 


14 


0.3393 


2.269 


4.33 


14 


0.3997 


1.927 


1.20 


15 


0.3253 


2.367 


0.91 


15 


0.3641 


2.115 


0.21 


16 


0.3194 


2.411 


0.07 


16 


0.3593 


2.143 


0.30 


17 


0.3116 


2.471 


0.41 










18 


0.3029 


2.542 


0.00 










19 


0.3146 


2.448 


0.01 



Table 13: Results of the quadratic fits to the p meson masses versus t . The inverse lattice 
spacing is defined by a -1 = 0.77GeV /m p (K c ). 



P = 5.85 


P = 6.0 


to 


m N > 


: 2 /df 


*o 


itln ; 


( 2 /df 


9 


0.6085 


3.20 


12 


0.4828 


0.03 


10 


0.6071 


0.37 


13 


0.4802 


0.24 


11 


0.6039 


0.83 


14 


0.4758 


0.79 


12 


0.5946 


2.15 


15 


0.4759 


1.88 


13 


0.5893 


0.37 


16 


0.4623 


0.05 


14 


0.5680 


2.28 


17 


0.4538 


2.00 


15 


0.5501 


1.59 


18 


0.4553 


2.28 


16 


0.5312 


0.49 


19 


0.4731 


0.10 


17 


0.5630 


2.55 


20 


0.4559 


0.10 



Table 14: Results of the linear fits to the nucleon masses versus t . 
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(3 = 5.< 


35 


(3 = 6.0 


to 


m A 




*o 


m A 


X 2 /df 


11 


0.6982 


0.03 


12 


0.6055 


0.01 


12 


0.6928 


0.01 


13 


0.6059 


4.10 


13 


0.6640 


0.05 


14 


0.6164 


1.78 


14 


0.6899 


4.02 


15 


0.6279 


0.50 


15 


0.5375 


34.96 


16 


0.6048 


0.29 


16 


0.4757 


23.96 


17 


0.6206 


0.12 








18 


0.6462 


0.17 








19 


0.5935 


0.34 



Table 15: Results of the linear fits to the A masses versus to- 



K 


Z v 


0.1440 


0.5121(9) 


0.1540 


0.5164(10) 


0.1585 


0.5126(36) 


0.1595 


0.5112(48) 


0.1605 


0.5101(76) 


Average 


0.5125(30) 



Table 16: Renormalization constants Zy for the local lattice current at (3 = 5.85 obtained 
in a previous work 0. 



(3 = 5.85 


(3 = 6.0 


K Fy T Fy P Fy C 


K Fy 1 ' Fy^ Fy C 


0.1440 0.2214 0.2600(33)1^ 0.1374 
0.1540 0.2211 0.2089(38)1? 7 0.1372 
0.1585 0.2094 0.1781(68)1^ 0.1299 
0.1595 0.1996 0.1658(84)1^ 0.1239 
0.1605 0.1885 0.1528(152)1^? 0.1170 


0.1450 0.1753 0. 1919(18)1^ 0.1210 
0.1520 0.1673 0.1558(17)1^ 0.1155 
0.1550 0.1544 0.1336(33)1^ 0.1065 
0.1555 0.1493 0.1276(43)1^ 6 0.1031 
0.1563 0.1382 0.1157(95)1^ 0.0953 



Table 17: p meson decay constants in lattice units. In parentheses are errors estimated by 
the jack-knife method. Errors given in the form lJ3er are f° r the fitting range dependent 
upper /lower bound. 
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(3 = 5.85 


P = 6.0 






pMC 






-pMC 


to 


Latt. 


Phys. 


X 2 /df 


Latt. 


Phys. 


X 2 /df 


to 


Latt. 


Phys. 


X 2 /df 


Latt. 


Phys. 


X 2 /df 


8 


0.169 


325 


53.7 


0.133 


257 


55.7 


11 


0.127 


297 


1.9 


0.107 


250 


2.0 


9 


0.157 


303 


6.0 


0.125 


o /in 

240 


7.0 


12 


f\ 1 O A 

0.124 


289 


1.7 


n i nr 

0.105 


O A A 

244 


1.9 


10 


0.156 


301 


11.1 


0.124 


239 


11.9 


13 


0.120 


280 


0.1 


0.102 


237 


0.2 


n 


0.145 


280 


5.8 


0.115 


222 


6.3 


14 


0.120 


278 


3.4 


0.101 


235 


3.6 


12 


0.141 


271 


0.0 


0.112 


216 


0.1 


15 


0.111 


259 


0.4 


0.094 


220 


0.4 


13 


0.130 


251 


2.3 


0.104 


201 


2.6 


16 


0.107 


248 


0.0 


0.091 


211 


0.0 


14 


0.148 


286 


0.3 


0.118 


226 


0.3 


17 


0.101 


234 


0.1 


0.085 


199 


0.1 


15 


0.116 


224 


0.0 


0.094 


180 


0.0 


18 


0.094 


220 


0.1 


0.080 


187 


0.1 


16 


0.105 


203 


1.1 


0.085 


164 


1.1 


19 


0.104 


243 


0.1 


0.088 


206 


0.1 



Table 18: Results of the linear fits to the p meson decay constants versus t , in lattice units 
and in physical units (MeV). 



P = 5.85 


P = 6.0 






= 12 


— t-^2 




^0 t m i n 


= 15 


t0 ty2 


K 




X 2 /df 






X 2 /df 


K 


rrii 


X 2 /df 


t x 2 


"I* 


X 2 /df 


0.1440 


1.0299(14) 


1.32 


8 


1.0304 


1.68 


0.1450 


0.8059(9) 


0.49 


12 


0.8068 


1.06 


0.1540 


0.6106(21) 


0.96 


7 


0.6117 


1.18 


0.1520 


0.4747(14) 


0.35 


7 


0.4767 


0.91 


0.1585 


0.3753(34) 


0.97 


5 


0.3774 


1.41 


0.1550 


0.2937(24) 


0.88 


6 


0.2967 


0.89 


0.1595 


0.3070(42) 


1.00 


4 


0.3097 


1.45 


0.1555 


0.2559(30) 


0.85 


6 


0.2593 


0.83 


0.1605 


0.2127(64) 


1.02 


4 


0.2175 


1.31 


0.1563 


0.1804(66) 


0.68 


5 


0.1897 


0.75 



Table 19: Pion masses determined from tx propagators. 



P = 5.85 


p = 6.0 


K fP s f'ps 


fPT fTP fMC 

J PS J PS J PS 


0.1440 0.1152 0.1443(25)1^ 
0.1540 0.0922 0.0929(22)11° 
0.1585 0.0732 0.0664(26)111 
0.1595 0.0677 0.0600(30)1^ 
0.1605 0.0597 0.0515(33)lf 


0.1450 0.0892 0.1030(10)1^ 0.0710 
0.1520 0.0713 0.0701(11)1J 8 0.0567 
0.1550 0.0566 0.0517(13)1$ 0.0450 
0.1555 0.0535 0.0482(15)1$ 0.0426 
0.1563 0.0462 0.0408(26)1^ 0.0368 



Table 20: Pseudo-scalar meson decay constants in lattice units. In parentheses are errors 
estimated by the jack-knife method. Errors given in the form Ij^r are f° r the fitting range 
dependent upper/lower bound. 
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Figure 1: Statistical errors in effective masses for the p meson at (3 = 6.0, K = 0.155 versus 
the bin size Nb- The errors are normalized by those for Nb = 1. 
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Figure 2: Statistical errors in effective masses for the p meson at (3 — 6.0, K = 0.155 versus 
the number of configurations N. The errors are normalized by those for N = N con f = 200. 
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Figure 3: Fitted mass mfu for the pion at (3 = 6.0, K = 0.155, obtained from one-mass fit 
to a range t — T/2 and the value of x 2 /df of the fit versus t. The error bars for mfu are 
statistical uncertainties estimated by the least mean square fit. Effective masses m e g with 
errors estimated by the jack-knife method are also given. 
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Figure 4: The same as Fig. |3| for the p meson. 
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Figure 5: The same as Fig. |3] for the nucleon. 
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Figure 6: Effective masses for the pion: a) f3 — 5.85, b) f3 — 6.0. The result of one-mass fit is 
reproduced by the solid line, dotted lines and dashed lines for the fitted mass, its statistical 
error and systematic upper/lower bounds, respectively. 
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Figure 7: The same as Fig. |6] for the p meson: a) (3 — 5.85, h) (3 = 6.0. 
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Figure 8: The same as Fig. ^|for the nucleoli: a) (3 — 5.85, b) (3 = 6.0. 
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Figure 9: The same as Fig. | for the A baryon: a) (3 = 5.85, b) j3 — 6.0. 
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Figure 10: Nucleoli to p mass ratio versus pion to p mass ratio squared. The errors shown 
are statistical only. The solid curve is obtained from phenomenological mass formulae . 
The dotted line is obtained by assuming that /m p and (m n /m p ) 2 are linear in the quark 
mass. The experimental value is marked with a star. 
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Figure 11: Masses of the ground state and the excited state for the p meson at (3 = 5.85, 
K = 0.1585 together with the value of x 2 /df °f the two-mass fits versus t m i n . The error 
bars are statistical uncertainties estimated by the least mean square fit. The result of the 
one-mass fit is reproduced by the solid line, dotted lines and dashed lines for the fitted mass, 
its statistical error and systematic upper/lower bounds, respectively. Note the difference in 
the scale of the plots for m and m 1 . 
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'lire 14: The same as Fig. |Tl| for the nucleoli at /3 = 6.0, K = 0.155. 
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Figure 15: Mass of the excited state of the p meson (denoted by p') at /3 = 5.85 versus 
1/K — 1/K C . The corresponding experimental values are marked with stars. The data for 
m p and m% are taken from the results of one- mass fits. With the scale of the plot, the results 
for nip from two-mass fits are indistinguishable from the one-mass fit results. 



1.5 



1.0 ■ 



0.5 ■ 



0.0 




0.00 0.05 0.10 0.15 0.20 
1/K- 1/K 



Figure 16: The same as Fig. [15] for (3 = 6.0. 
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Figure 17: Mass of excited state of the nucleoli (denoted by N') at (3 — 5.85 versus 
1/K — l/K c . The experimental values for masses of the nucleon (N) and its excited state (TV 7 ) 
are marked with stars. The data for m^, m p and are taken from the results of one- mass 
fits. With the scale of the plot, the results for from two- mass fits are indistinguishable 
from the one-mass fit results. 
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Figure 18: The same as Fig. 17 for (3 = 6.0 
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Figure 19: Linear extrapolations of hadron masses at (3 = 5.85 to the chiral limit. The open 
circles at zero quark mass are extrapolated values. The errors shown are statistical only, and 
do not include the systematic errors discussed in the text. 
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Figure 20: The same as Fig. 19 for (3 = 6.0. 
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Figure 21: p meson masses at (5 — 5.85 obtained from the fit with various to together with 
linear extrapolations of these data. The open circles are extrapolated values. The errors 
shown are those estimated by the least mean square fits. 
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Figure 22: The same as Fig. but for /3 = 6.0. 
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Figure 23: a 1 determined from linear fits to all possible combinations of the p masses 
obtained by varying t from t x 2 to 18. 
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Figure 24: Ratio of the vector meson decay constant to the vector meson mass, for the three 
choices of renormalization constants discussed in the text. The errors shown are statistical 
only and are estimated by the jack-knife method. The corresponding experimental values for 
vector mesons are marked with stars. The value of mps for the strange quark is estimated 
by phenomenological mass formulae |15| using my = = 1019 MeV. 
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Figure 25: Ratios of the (f) meson decay constant and the K meson decay constant to the p 
meson mass versus the p meson mass in lattice units. The errors in our data are statistical 



only. The GF11 data are taken from ref. [21]. The corresponding experimental values are 
marked with stars. 
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Figure 26: Linear extrapolations of vector meson decay constants, for the two choices of 
renormalization constants discussed in the text. The open symbols at zero quark mass are 
extrapolated values for (3 = 6.0. The errors shown are statistical only. The experimental 
value for the p meson is marked with a star. 
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Figure 27: The same as Fig. 25 for the p meson and the pion. 
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Figure 28: Ratio of the pseudo-scalar meson decay constant to the vector meson mass, for 
the three choices of renormalization constant discussed in the text. The errors shown are 
statistical only and are estimated by the jack-knife method. The corresponding experimental 
values for pseudo-scalar mesons are marked with stars. 
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Figure 29: Linear extrapolations of the pseudo-scalar meson decay constant, for the two 
choices of renormalization constant discussed in the text. The open symbols at zero quark 
mass are extrapolated values for (5 = 6.0. The errors shown are statistical only. The 
experimental value for the pion is marked with a star. 
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